A resolution of the transition to turbulence paradox 

Understanding the transition to turbulence has been impeded by a modehng oversight 
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Despite being around for over a century, the tran- 
sition to turbulence problem remains central in fluid 
dynamics. This phenomenon was apparently known 
to Leonardo da Vinci [15J, who in 1507 introduced 
_ the term "la turbulenza", and nowadays it has an 
impact on practically every field ranging from astro- 
physics and atmosphere dynamics to nuclear reactors 
and oil pipelines. Beginning with the systematic ex- 
^ perimental studies in a pipe by Osborne Reynolds 
[l] in the 1880s, it is known that the flow becomes 
turbulent at finite flow rate, usually measured by 
Reynolds number Re — LU /v (see the definitions in 
^ figure [1]). Similar observations have been made in 
K^> other flows, in particular Couette flow the flow be- 
I tween two plates moving parallel to each other, cf. flg- 
jSure 
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1(a) , where the transition is observed at Re : 
^ (see, for instance, [2]). Reconciliation of these experi- 
O inental observations with theory [3l I14| failed because 
^ the eigenvalue analysis of the linearized Navier-Stokes 
(~| equations (NSEs), which govern fluid motion, yields 
I ^E igenvalues A in the left half-plane at all values of 
Rc, which implies that all small initial disturbances 
<^ should decay exponentially like e'^*, as time increases 
and thus one should have stability. The basic math- 
— . ematical setup in classical works [3] treats the stabil- 
ity problem in an infinite channel x e ( oo,+cxd), as 



C ■ in figure 1(a) , In this work we demonstrate that this 



QQ infinite channel assumption is a sticking point that 
has prevented one from understanding the primary 
^ instability in the transition to turbulence. Our anal- 
ysis on semi-infinite channel-domain x G [0, +oo), cf. 



*^ figure 1(b) , which is more relevant to the way exper- 



iments are usually done, predicts instability and thus 
explains many important features of these phenom- 
ena in a simple and basic way. 

First, let us recall the NSEs, which can be written as an 
evolution equation for the velocity field u. 



du 

— = ^u 

dt 



N(u), 

with the linear and nonlinear terms given by 

Au = P [-U ■ Vu - u • VU + Re-^Au] 
N{u) = -P [V • (u u)] , 



(1) 

(2a) 
(2b) 



where P is the projection on the space of divergence free vector 
fields and U is the base flow. In the class of weak solutions, 
Romanov |14| established linear and nonlinear stability, based 
on the absence of linearly unstable eigcnmodes. of the Couette 
base flow U = {y,0) on an infinite domain x G (— cxd, +cxd). 
This, of course, leads to a conceptual difficulty since instabil- 
ity is observed in experiments. 
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(a) Infinite channel setup, x £ 
(— oo, +00). 
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(b) Semi-infinite channel setup, x G 
[0,-1-00). 

Figure 1: Couette flow of a fluid with kinematic viscosity v 
in a channel of width L; the base flow is U = (y, 0). 

The failure of the hydrodynamic stability theory based 
on ([1]) to predict the transition to turbulence motivated vari- 
ous alternative explanations, including the idea of a very small 
basin of attraction of the stable base flow and the transient 
growth idea [H [5] . It is also understood that the transition 
to turbulence belongs to a general class of counter-intuitive 
dissipation-induced instabilities based on the recent theory 
[TOl [TTJ [12] . While all these approaches are still in develop- 
ment, it is worth mentioning the line of logic of the transient 
growth concept. Namely, based on the ansatz that the non- 
linear terms (j2bp of the NSEs ([T|) arc energy conserving and 
since the linear terms (|2ap can produce energy only transiently 
in time, then the transient growth is the only explanation of 
the fact that we observe non-zero deviations from the lami- 
nar base state U in the aforementioned flows. The transient 
growth itself is related to the sensitive and non-normal nature 
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of (Hil) g], i.e. AA* ^ A* A, where A* is the operator adjoint 
to A. While this picture has been appended with various dy- 
namical systems scenarios [5], such as chaotic saddles, and it 
also appears to be useful as a transient effect , there is still 
no theory which would be able to predict a transition robustly. 
Here we propose a more direct resolution of this long-standing 
problem by demonstrating the existence of linearly unstable 
eigenmodes. 

Before giving a resolution of the mismatch of theory and 
experiment, we would like to remind the reader that the fact 
of instability implies an existence of at least one eigenmode 
fx{x) such that the corresponding eigenvalue A has a positive 
real part. Perhaps because of the translational invariancc of 
the base state U ~ (y, 0), cf. figure 1(a) or for convenience. 



domain is versus —ik, where k is the wavenumber. 



the stability of Couette flow has always been studied on an 
infinite domain, — oo < x < +00. However, if one recalls 
the way the experiments on the transition are usually done, 
i.e. one introduces disturbances at some fixed inlet location, 
say a; = 0, and observes how they evolve downstream, then 
it becomes clear that the semi-infinite domain, x G [0, -l-oo), 



as in figure 1(b) is more relevant as a mathematical ideal- 
ization. Then, it is convenient to study the linear eigenvalue 
problem - the classical Orr-Sommerfeld (OS) equation - by 
assuming that the disturbance eigenfunction is of the form 
~ axtj.{y)e^*e~^^^, where A G C is the eigenvalue and /i > is 
an analog of the wavenumber: 



^ + /i^ ) a{y) = 0, (3a) 



-1,1 



(3b) 



where we dropped the indeces A and fi. The basis functions 
e~^^, fi>0, arc clearly not members of the space of bounded 
functions on the whole real line, x G (—00, +00). which are 
used in the classical analysis of this problem, but they do be- 
long to the space of functions bounded on the semi-infinite do- 
main, X G [0, +cxd). Therefore, these eigenfunctions were not 
captured in the traditional approach. Alternatively, equation 
^ could be treated with cosine and sine Fourier-transforms, 
which would lead to the same results but with a considerably 
more complicated version of ([3]). We choose to work with 
exponential functions e^^^', since our goal here is to demon- 
strate the existence of unstable eigenmodes in a direct way; 
other reasons will be clear from the subsequent discussion. 
Also, this choice of functions leads to the OS equation anal- 
ogous to the classical one [5]: the only difference of ^ from 
the classical case of the OS equation ^ studied on an infinite 
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(a) max_Re(A); o: Rc = lO"*, X: Re = 10^, 
□: Re = 10^, +: Re = 10^. 




(b) Critical Reynolds number Rcc(/j). 
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(c) Asymptotics max/?e(A) for Re = 10*. 
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(d) Spectrum: distribution of eigenval- 
ues An, ?i G N of l(3]l in the complex 
plane (eigenvalues continue to the neg- 
ative part of the real axis) for = 1 and 
Re = 70 (red dots), Re = 1000 (green 
dots), and Re = 5000 (blue dots). 
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Figure 2: Stability picture of the Couettc flow. 

The eigenvalue problem ([3]) is solved numerically by ex- 
panding a{y) in Chebyshev polynomials [5]. As one learns 
from figure m the stability picture on a semi-infinite domain 
is in sharp contrast to what one has on an infinite domain [3] , 
but conforms well with the usual intuitive understanding of 
instability phenomena: for some values of ^ there are eigen- 
values with positive real part and thus Couette flow is abso- 
lutely unstable. In fact, since in general all values of /i may 



be present in a real flow, then figure 2(a) suggests that the 



transition in the Couette flow is not a critical phenomenon. 



Indeed, figure 2(c) indicates that the instability in the Cou- 



ettc fiow is in fact a short-wave instability, since the value of 
maxi?e(A) is increasing with /i, and larger /i means that the 
disturbance is localized around the inlet. It is also notable 
that the structure of the eigenfunctions corresponding to the 
leading eigenvalues increases in complexity with increasing /i, 
as illustrated in figure [3l which may explain the tangled flow 
picture observed experimentally in the Couettc flow: see |16| 
and references therein. However, if in a particular experiment, 
the admissible magnitudes of ^ are restricted to a range of 
small values, then one can observe critical phenomena, as in 
figure [2(b)| These are clearly of Hopf bifurcation type, com- 
mon in various fiuid dynamics problems |13] . as follows from 
the distribution of eigenvalues in the complex plane in figure 



2(d) which illustrates that the leading (rightmost) eigenval- 



ues cross the imaginary axis as Reynolds number increases. 
In this case one can expect that the leading eigensolutions are 
the usual ToUmien-Schlichting waves [l0l[3] appearing via the 
Hopf bifurcation. We have to stress, however, that in general, 
if all values of fi are present, then the transition to turbulence 
is not a critical phenomenon and thus not a Hopf bifurcation, 
similar to the Rayleigh- Taylor instability, i.e. instability of a 
heavy fluid accelerating into a light one. 




Figure 3: The eigenfunctions corresponding to the rightmost 
eigenvalue for the Couette flow at Re ~ 5000: red curve corre- 
sponds to ^ = 1 and Rc Amax — 38.53, blue curve corresponds 



to ^ = 40 and Re Amax — 0.76. 

There is much more to this stability picture and a lot 
remains to be understood about the properties of equation 
Q, as well as a full function-analytic nonlinear analysis of 
([T]) with careful treatment of the inlet boundary conditions 
is needed similar to [12] . However, the key feature - the ex- 
istence of unstable eigenmodes - originating from the semi- 
inflnite domain setup is well illustrated above. Analogous 
computations performed for the plane Poiseuille flow also re- 
vealed a similar instability picture, which suggests that the 
right mathematical setup used above is a universal expla- 
nation of the transition to turbulence in the aforementioned 
"troublesome" flows. 

This counter-intuitive difference in the stability results 
between the semi-infinite and infinite domains can be appre- 
ciated with the assistance of the sketch in figure [D 

fix) 
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Figure 4: Eigenmodes in the stability analysis: semi- versus 
infinite domain. 

Namely, if, for example, one restricts (eigen-) functions 
f{x) to be bounded for all x as dictated by the fact that 
the physical solution should be bounded, then eigenfunctions 
defined on x G (— cx), +oo) are more restrictive compared to 
eigenfunctions deflned on x G [0, +oo). Indeed, if one can 
construct a function /+ bounded on a; G [0, -f-oo) which also 
satisfies the OS equation (|3ap . then continuation of this func- 
tion onto x G (— cxo, 0] may lead to an unbounded function /_, 
as dictated by the structure of the linear operator (|2ap and as 



illustrated in figured Our exponential eigenfunctions e^*^^ 
are a good example of functions bounded on the right half-line 
and unbounded on the left half-line, while still satisfying the 
OS equation ([3]). These observations can be illustrated with 
the following elementary eigenvalue problem: 

d^<j>{x) d^<f>{x) 



X- 



0, 



(4) 



dx^ dx^ 

which clearly contrasts the problems on infinite and semi- 
infinite domains: 

• a; G (— oo,+C)o) and (/>(a;) G L2'. applying Fourier trans- 
form we get A = — ifc, fc G M, i.e. marginal stability. 

• a; G [0, -l-oo) and \4>{x)\ < oo: instability is present since 
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there are eigenfunctions 
values \ = iJ,. 



/i > 0, with eigen- disappear, since disturbances are non-zero at the inlet and, if 
they lead to an instability, do not necessarily decay at infinity. 



Note that the above arguments also explain the sensitiv- 
ity of the experimentally observed critical Reynolds number 
Rcc to the properties of disturbances at the domain inlet, 
X = 0: while their amplitudes do not play a role in view of 
the linearity of the problem, gradient-like properties of the 
disturbances do! Namely, by varying gradient-like proper- 
ties of disturbances at x = effectively changes the bound- 
ary conditions at .t = and thus the size of the eigenspace. 
Since restricting the domain to x £ [0, -f-oo) enlarges the 
function space, one can expect that the spectrum enlarges 
as well and may lead to instabilities. This simple fact, as 
we saw above, explains the mechanism behind the transition 
to turbulence! While the above analysis has been conducted 
on a semi-infinite domain with open inlet and outlet bound- 
aries, which demonstrated the existence of absolutely unstable 
eigenmodes, the latter apparently also exist on finite length 
domains with open inlet and outlet boundaries. 

Finally, note that there is nothing unusual about increas- 



ing growth rate for increasing /.t, observed in figure 2(c) which 



is a common feature in many fundamental hydrodynamic in- 
stabilities in the short-wave limit, such as the Rayleigh- Taylor 
instability. Of course, in the nonlinear setting one does not 
observe infinite growth rates, because they are suppressed by 
the nonlinear terms, which play a stabilizing role opposing the 
formation of singularities even in the case when they are en- 
ergy conserving [3]; in our case the nonlinear terms can also 
be dissipative. The latter fact is in apparent contrast with 
the above mentioned ansatz of the transient growth story, i.e. 
that the nonlinear terms are energy-conserving. This ansatz 
is valid if the disturbance field in the Couette flow problem 
is considered on an infinite domain with the boundary con- 
ditions at X = ±oo corresponding to decay to zero, but not 
on a semi-infinite domain. Indeed, multiplying with 
and integrating over the semi-infinite strip-like flow domain 
we arrive at the Reynold-Orr equation for the kinetic energy 
E{t) = ||u||V2: 



— = / TliUipds - 

an 



(5) 



-I- {Dij,UiUj) + 



an 



ds 



+ n / rijUjUiUi dx, 
2 Jo 



where n is the outward normal of d^l. If the domain H, is 
unbounded and open, as in the Couette or pipe flows, then 
the effect of the nonlinear terms (cubic term in ([5])) does not 
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